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Abstract 

We experimentally studied the effect of positional disorder on a Josephson junction array with 
/ = n, | +71, or | + n flux quanta per unit cell for integral n. This system provides an experimental 
realization of a two-dimensional XY model with random phase shifts. Contrary to many earlier 
numerical and analytical investigations, our results suggest that low-temperature superconductivity 
is never destroyed by positional disorder. As the disorder strength increased, the Kosterltz-Thouless 
(KT) type order in the / = and ^ systems changed to a non-KT type order with a long-range 
phase coherence, which persisted even in the maximal disorder limit. A possible finite-temperature 
glass transition is discussed. 

PACS numbers: 64.60.Cn, 74.81.Fa, 74.25. Qt, 75.50.Lk 
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The XY model with quenched random phase shifts, also known as the random-gauge 
XY model, has been intensively studied with a focus on the maximal disorder limit as a 
■node, for a vo rteX g.ass phase ft of h lgh -te m pe r atu rc superconductors in a magnet.c Md. 
Despite a great deal of effort for the last two decades, the model in two dimensions (2D) 
remains not well understood. Conflicting evidences, mostly numerical or analytical ones, 
forbid a consensus about the possibility of a finite-temperature glass transition in 2D. The 
purpose of this work is to see experimentally how the low-temperature ordered phase of the 
2D random-gauge XY system evolves with the disorder strength increased. 

The random-gauge XY model is described by the Hamiltonian 

H = -JY,cm(9 i -9 j -A ij ), (1) 

where 9i is the phase of the order parameter on site % and the quenched random phase shifts 
A^s are uncorrelated and distributed uniformly in [A°a — m, A°a + rir] with < r < 1. The 
case with maximal disorder corresponds to the disorder strength r = 1. The effect of random 
phase shifts in 2D was investigated primarily with respect to the Kosterlitz-Thouless (KT) 
transition for / = 0, where frustration / = (l/27r) The directed sum of the is 

around a unit cell. A prevailing conclusion in many numerical and analytical studies [2j for 
/ = is that there exists a KT transition at finite temperatures when r < r c « 0.37 and 
the case with maximal disorder, termed the gauge glass, has no ordered phase at any finite 
temperature. The lack of ordering at finite temperatures has been proven rigorously in 
the limit of maximum disorder. Nevertheless, numerical evidences against this conclusion 
have been also provided by several different research groups . The counter evidences 
point to a finite-temperature superconducting glass or de r in the gauge glass limit. It has 
been shown p| numerically that a breaking of ergodicity due to the large energy barrier 
against vortex motion may open up the possibility of a finite-temperature glass transition 
which was denied by the analytical equilibrium studies (^J. The confused situation is much 
the same for / = |, where the XY system without disorder has the Ising-like structural 
order as well as the KT-like phase coherence at low temperatures |jj . Along with arguments 
[s| for the persistence of both the Ising-like order and the KT-like order for at least weak 
disorder strength, it has been also suggested |9 ; ] that the critical disorder for / = ~ is zero, 
or that any amount of disorder, no matter how weak it is, destroys the Ising-like order and 
the phase coherence even at zero temperature. 



The disagreement among the analytical and numerical results may be settled by exper- 
imental evidences for the random-gauge effect. The random-gauge XY model in 2D can 
be most closely realized by a Josephson junction array (JJA) with positional disorder in a 
uniform magnetic field [10]. For a JJA, J is the Josephson coupling energy per junction and 
given by J = hi c /2e with i c the single-junction critical current. 9i is the phase of the super- 
conducting island at site i and Aij = (27r/0o) H A • dl is the integral of the vector potential 
across a junction with 0o being the flux quantum. The frustration / is equal to the mean 
number of flux quanta penetrating each cell. Random phase shifts can be introduced by 
randomly displacing the centers of superconducting islands while keeping the junctions on 
their regular positions in a periodic lattice. For a square array with frustration / and with 
the random displacements in x and y directions distributed uniformly in the range [—A, A] 



with A in units of the 
relation r m 2.2 / A JlO 



attice constant, the corresponding disorder strength is given by the 



111 ]. Since the XY Hamiltonian is invariant under the transformation 
of / to f + n for an integer n, the disorder strength of a JJA for a specific / can be controlled 
by changing either n or A. In this paper, we present an experimental investigation of a posi- 
tionally disordered JJA with f = n, | + Ji, or | + n for various n's. Transport measurements 
revealed that the three different systems superconduct at sufficiently low temperatures for 
any disorder strength r and, despite the essential differences among them at weak disorder, 
belong to the same universality class in the limit of maximal disorder, proposing the possi- 
ble occurrence of a finite-temperature glass transition. For / = and |, the nature of the 
superconducting transition was found to change from a KT type to a non-KT type with r 
increased. 

Experiments were performed on a square array of 200x800 Nb/Cu/Nb Josephson junc- 
tions with a controlled amount positional disorder, a section of which is shown in Fig. ^a). 
Nb islands were disposed on a 0.3-/iin-thick Cu film periodically with a lattice constant 
of 13.7 jttm, a junction width of 4 /xm, and a junction separation of 1.4 /jm. The random 
displacements of the centers of the crosses in x and y directions were uniformly distributed 
in the range [—A, A] with A = 0.20 in units of the lattice constant. For an ideal JJA 
with zero junction width, this amount of displacement corresponds to the disorder strength 
r rs 0.44/. However, the actual disorder strength for our sample with a finite junction 
width must be weaker than designed, due to the asymmetric diamagnetic responses of the 
distorted superconducting islands which inhomogeneously changes the effective area of each 
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cell to reduce the amount of positional disorder |12j. If the relocation of the effective center 
of the Nb island is limited within the middle square of the cross, the smallest possible value 
of the effective displacement Ag/j for A = 0.20 is ~ 0.05. The actual amount of A e yy of 
the sample will be determined in the stage of comparing the experimental results with the 
numerically proposed phase diagram. In the absence of a magnetic field, i c and J (= hi c /2e) 
can be obtained by extrapolating the low-temperature i c vs T data, obtained from the / 
vs dV/dl curves, by using de Gennes formula for proximity-coupled junctions in the 
dirty limit. In the presence of a magnetic field, i c depends on the magnetic flux through 
the junction i c {B) = i c (0) sin (7r0/0o ) /(y </>/0o )• Since it is not possible to directly 

measure i c {B) of a positionally disordered array and i c (B)/i c (0) (= sin (7r0/0 o )/ (^0/00)) 
depends mainly on the geometry of the junction, i c (B) of the sample was determined by 
combining i c (0) of the sample and i c (B)/i c (0) (= sin(7r0/0o)/(7T0/0o)) of another array with 
the same junction dimensions and properties but no positional disorder. The resistance (R) 
and the current- volt age {IV) characteristics were measured by employing the phase- sensitive 
voltage-signal-detection method using a lock-in voltmeter and a square-wave current at 23 
Hz. The magnetic field or the frustration / was adjusted from the / vs R curve of the sample 
exhibiting distinct resistance minima at integral /'s. More details of the experiments are 
described in Ref. 

Figure ^b) shows the resistive transition of the sample for three different sets of frustra- 
tions / = n, \ + n, and \ + n. Although the transition broadens appreciably with increasing 
n, it is evident that for all the frustrations, the transition to a zero-resistance state takes 
a place at a finite temperature. Figure El illustrates the / dependence of the supercon- 
ducting transition temperature T c determined at R = 5 x 10 -6 Q. Note that T c is plotted 
in units of J/ks in order to scale out the temperature and magnetic-field dependences of 
the critical current i c (T, B) and to make a direct comparison of the T c with the numerical 
results possible. The T c 's determined from the resistance measurements agree within error 
bars with those independently determined from the analysis of the IV characteristics (to 
be commented on in detail below). We first consider the / = case. As n (or the disorder 
strength) increases, T c decreases at first from 0.84J/ks and then stays at ~ 0.35 J/ ks for 
f > f c ~ 3. The phase line never ends at T = 0, contrary to many earlier investigations 
P, [J. A similar phase boundary with f c ~ 5 was observed for the sample with A = 0.15. 
The nature of the superconducting transition could be identified from the development of 
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the IV characteristics with temperature. Measurements of the IV characteristics were car- 
ried out for six different frustrations / = 0, 2, 3, 4, 6, and 8. Shown in Fig. El are those 
for f — 2 and 6. The IV curves were obtained by averaging 15-240 measurements for each 
current. For / = 2, the IV curves below the transition display a power-law behavior (a 
straight line in a log/ vs \ogV plot or a flatted curve in a log/ vs /(log V)/d(\ogI) plot) at 
low currents, indicative of the KT-type order with a quasi-long-range phase coherence and 
vortex-antivortex pairs as dominant excitations. The fast drop of the slope d(log V)/<i(log/) 
of the IV curves with temperature above the transition is also characteristic of the KT tran- 
sition. The IV curves for / = 6, on the other hand, have an exponential form (a convex 
upward curve in a log-log plot) below the transition, indicating a non-KT type order with 
a long-range phase coherence. As shown at the bottom of Fig. the TV curves for / = 6 
can be fitted to the Fisher-Fisher-Huse scaling form Q for a non-KT type superconducting 
transition in 2D, TV j!\T — T c \ zu = e +(I/ T\T — T c \ u ), where z is the dynamic exponent and 
v is the correlation length exponent [15j. The critical exponents from the scaling analysis 
are v = 2.0 ± 0.3 and z — 1.8 ± 0.3 for / > 3, insensitive to /. Crossing over from the 
KT-like behavior to the non-KT-like one arose at / ~ f c (~ 3), the frustration above which 
T c became approximately independent of /. The evolution of the IV characteristics with 
/ discloses that the nature of the superconducting transition for / = changes from a KT 
type at weak disorder to a non-KT type at strong disorder. These experimental findings 
except the values of T c and v at strong disorder are in a good agreement with the numerical 
results reported in Ref. |5|. Comparing the phase boundary for / = in Fig. El with the 
numerically proposed phase diagram in Ref. p, we find that f c ~ 3 for our sample corre- 
sponds to the critical disorder strength r c « 0.37 and thus / ~ 8 to r = 1. This implies a 
finite-temperature superconducting transition even in the gauge glass limit of r = 1. The 
correspondence of f c ~ 3 with r c « 0.37 also gives ~ 0.06 for A e ff of the sample. 

The effect of positional disorder on the / = \ system which has the KT-like phase 
coherence and the Ising-like structural order in the absence of disorder were found quite 
similar to the observed for / = 0. As shown in Fig. with increasing n, T c of the / = | 
system drops fast from 0.45// fee (the T c in the absence of disorder) to ~ 0.25 J / k B and then 
increases slowly until it reaches ~ 0.3J/k B . The IV characteristics for / = ~, 1^, 3|, 5|, 
and 7|, some of which are shown in figures 0] and Efb), indicate that similarly to the / = 
case, the nature of the superconducting transition changes from a KT type at weak disorder 
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to a non-KT type at strong disorder. The critical exponents for the non-KT type transitions 
are v = 2.0 ± 0.3 and z — 1.8 ± 0.3. The major difference from the / = case is that the 
crossing over appears at lower f c < 1.5. 

Dissimilarly from the / = and | systems, the / = | system without disorder has a 
superconducting vortex-solid phase with a long-range phase coherence at low temperatures 
and experiences a melting transition driven by domain wall excitations [16]. The IV char- 
acteristics for / = |, l|, 3|, 5|, and 7|, some of which are shown in figures |U and Efb), 
reveal that for / = |, the introduction of positional disorder does not alter the essential 
features of the low-temperature ordered state and phase transition. The phase diagram of 
Fig. |21 displays T c for / = | as a monotonously increasing function of the disorder strength. 
The scaling analyses of the IV data present 1.8-2.0 for v and 1-1.8 for z, dependent on the 
disorder strength. 

The experimental results for the three different sets of frustrations demonstrate that as 
the disorder strength (or the amount of the random phase shifts) increases, a KT type 
order of a 2D XY system weakens and is eventually replaced by a non-KT type order 
with a long-range phase coherence, which never weakens even at the maximum disorder 
strength. A remaining important question to be answered through this work is whether 
the finite-temperature superconducting phase in the maximal disorder limit is a glass in 
nature. Figure Efa) exhibits the variation of the sample resistance with / at T = 4.20 K. 
The resistance oscillation fades away as / approaches the maximal disorder limit (/ ~ 8). 
Three sets of IV curves for / = 8, 7|, and 7§ plotted in a single figure [Fig. EJb)] also shows 
that the IV characteristics near maximal disorder are independent of frustration, excluding 
the slightly /-dependent T c . Consequently, the three sets of IV data collapse onto the same 
curves with / and V scaled by the Fisher-Fisher-Huse scaling form with the same v (= 2.0) 
and z (= 1.85). This suggests that the 2D XY systems with the maximal disorder strength 
belong to the same universality class, regardless of the amount of frustration imposed on the 
system, or that the natures of the superconducting state and phase transition of a 2D XY 
system become independent of / in the gauge glass limit, despite the essential differences 
among the systems with different /'s at weak disorder. Although the data do not provide a 
full account for the nature of the low-temperature superconducting state, these features of 
the data seem to propose the possible occurrence of a finite-temperature glass transition in 
the gauge glass limit. 
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FIG. 1: (a) Photograph of the positionally disordered array with A = 0.20. (b) Resistive transition 
of the disordered array for / = n, \ + n, and | + n. 

FIG. 2: Superconducting transition temperature T c in units of J/ks for three different sets of 
frustrations f = n (filled circles), \ + n (open squares), and § + n (filled triangles), determined 
from the resistance measurements. 

FIG. 3: / vs V curves in a log-log scale for / = 2 at T = 4.200 to 5.200 K and for / = 6 at T = 
3.800 to 4.900 K; the slope d(log V)/d(log I) of the / vs V curves as a function of /; and scaling 
plots of the I vsV curves at 12-18 different temperatures. The dashed lines in the I vs V plots are 
drawn to show where the phase transition occurs. The insets of the scaling plots show the values 
of T c , v, and z used to scale the data. 

FIG. 4: Development of / vs V curves with temperature for / = |, 3|, |, and 3|. The dashed 
lines are drawn to show where the phase transition occurs. 

FIG. 5: (a) Frustration dependence of the resistance of the sample with an excitation current of 30 
fiA at T = 4.20 K. (b) I vs V curves for three different frustrations / = 8, 7|, and 7| in a single 
plot. 
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